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Abstract 



O |. We theoretically investigate the quantum noise properties of the pulse pumped high gain fiber 

rj ' optical parametric amplifiers (FOPA) by using the Bogoliubov transformation in multi-frequency 

^ ' modes to describe the evolution of the non-degenerate signal and idler twin beams. The results 

show that the noise figure of the FOPA is generally greater than the 3 dB quantum limit unless 

(N . 

^ \ the joint spectral function is factorable and the spectrum of the input signal well matches the gain 

OO ' 

spectrum in the signal band. However, the intensity difference noise of the twin beams, which 

m ■ 

lO ■ weakly depends on the joint spectral function, can be significantly less than the shot-noise limit 

f"^ \ when the temporal modes of the pump and the input signal are properly matched. Moreover, to 

cn ■ 

closely resemble the real experimental condition, the quantum noise of twin beams generated from a 
broadband FOPA is numerically studied by taking the various kinds of experimental imperfections 

X. 

^ . into account. Our study is not only useful for developing a compact fiber source of twin beams, 

but also helpful for understanding the quantum noise limit of a pulse pumped FOPA in the fiber 
communication system. 
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I. INTRODUCTION: 



It is well known that the fiber optical parametric amplifiers (FOPAs), employing x^^^ 
nonlinearity based four wave mixing (FWM) to transfer energy from one or two strong pump 
fields to a weak signal field, can be used in the conventional optical communication systems 
for signal-processing applications, such as optical amplification, phase conjugation, optical 
limiters, 3R generators, and time-domain de-multiplexer etc. \Mi^- ^^ fact, similar to the 
crystals based optical parametric amplifiers (OPAs) [sHsl , FOPAs are also candidates 



X 



(2) 



for generating non-classical light for quantum communication |9|-|ll|. 

In quantum communication, information can be encoded discretely on each photon or 
continuously in quadrature-phase amplitudes of an optical field. The former requires single- 
photon sources that can be generated with weak nonlinear interaction [12], but the latter 
relies on strong interaction to generate quantum correlation between the amplitudes of opti- 



cal fields for entanglement 



a 



131] . Moreover, continuous variable (CV) entanglement offers 



several advantages over its discrete variable counterpart. The main advantage is the ex- 
perimental feasibility in unconditional production of entangled states which in turn enables 
quantum information protocols to be accomplished unconditionally 13]. 
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14J-|l8J. In order 



The past decade has seen the growing interest in developing nonclassica 
the FWM in fibers, because they are compatible with the fiber network 
to enhance the FWM and to suppress the Raman scattering [l9[, a pulse pumped FOPA 
is preferred because the peak power of pulsed laser is very high even at a modest average 
power. However, for the majority of the experiments performed to date, the FWM in fiber is 
in the low gain regime, and the fiber based nonclassical light sources with vacuum injection 
are in the domain of a few photons and thus can be only employed for discrete variable 
information encoding. In order to develop the fiber based source of the CV nonclassical 
light, the pulse pumped FOPA needs to be operated in the high gain regime. 

Up to now, only a couple of experiments reported the generation of CV nonclassical 



light via the pulse pumped FWM in fibers [9|, |20|. In 2001, Sharping and co- works realized 
the pulsed twin beams by using the FOPA for the first time |9|. The amplified signal 
and generated idler beams bear a strong quantum correlation in intensity (photon-number) 
fiuctuations in the sense that the quantum noise level of their intensity difference is lower 



than the shot-noise hmit (SNL). In that experiment, the gain of FWM is less than 3 and 
the observed intensity difference noise of twin beams is only about 1.1 dB (2.6 dB after 
correction for losses) lower than the SNL. Recently, using the FOPA with photon number 
gain of about 16 dB, our group demonstrated the twin beams with the intensity-difference 
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The results indicate 



22| . the pulse pumped high gain 



noise below the SNL by 3.1 dB (10.4 dB after correction 

that compared with its x*-^"* crystal counterparts [Tl, |2l|, 

FOPA is a simple and promising system for generating the CV nonclassical light, which 

has the potential applications in realizing the quantum communication protocols and in 

performing the high precision measurement. 

There have been a few theoretical papers analyzing the generation of pulsed CV nonclas- 
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24| . For the pulse pumped high 



sical light via the high gain x*^^^ crystals based OPAs 
gain FOPAs, however, the investigation is mostly within the domain of classical character- 
istics |2|, y] and the quantum theory has not been worked out yet. This is different from 
the case of continuous wave (CW) or quasi-CW pumped high gain FOPAs |2|, [lO|, 111]. In 
order to understand the high gain FOPA for further reducing the intensity difference noise 
of pulsed twin beams, we focus on developing a multi-mode quantum theory in this paper. 
The rest of the paper is organized as follows. After briefly introducing the quantum 
description of a single-mode FOPA in Sec. II, we deduce a multi-mode quantum theory of 
a pulse pumped non-degenerate FOPA in Sec. III. The evolution of non-degenerate signal 
and idler beams is described by the Bogoliubov transformation in a multi-frequency mode 
form. According to the Hamiltonian of FWM in fibers, we derive the four transformation 
functions of the Bogoliubov transformation in the analytical form of infinite series. The 
calculation shows that the noise performance of the high gain FOPA highly depends on the 
joint spectral function (JSF). In Sec. IV, we show that the FOPA with factorized JSF can 
be described by a quantum model of single-mode FOPA, which is similar to the model given 
in Sec. II, because both signal and idler beams can be in single temporal mode, respectively. 
The results illustrate the validity of the multi-mode theory. In Sec. V, we focus on analyzing 
the FOPA with a very broad gain bandwidth. In this situation, the JSF is non-factorable. 
The calculated results show that although the pulse pumped high gain FOPA generally has 
a noise figure greater than the 3 dB quantum limit, the intensity difference noise of the twin 
beams, weakly depends on the JSF, can be significantly reduced to lower than the SNL. To 
closely resemble the real experimental situation, the dependence of the intensity difference 



noise of twin beams is numerically studied by considering the influences of experimental 
imperfections, including the quantum efficiencies of detectors, the collection efficiency of 
twin beams, and the excess noise of input signal. Finally, we conclude in Sec. |VIl 

II. QUANTUM DESCRIPTION OF A SINGLE MODE FOPA 



We start with the quantum description of single-mode linear FOPAs, whose pumps are 
single frequency lasers with negligible depletions. In this section, we will introduce the 
Bogoliubov transformation in the single frequency mode form and the analytical expressions 
of key parameters of the FOPA, including the photon number gain, the noise figure, and the 
intensity difference noise of the twin beams. 

Figure [Il^a) is a typical diagram of the non-degenerate phase-insensitive FOPA, which 
is referred to as the FOPA hereinafter for brevity. When the wavelengths of the pump 
are properly set, the phase matching of FWM in fiber is satisfied. In the FWM process, 
two pump photons at frequencies Upi and ujp2 are simultaneously scattered into a pair of 
signal and idler photons at frequencies Us and Wj via the x^^^ nonlinearity, and the energy 
conservation relation Upi + ujp2 = Ug + uJi is satisfied. Since each generated signal photon is 
always accompanied by the birth of an idler photon, the two form correlated twin beams. A 
weak input is injected into the FOPA from either the signal or the idler field. For convenience, 
we refer the field with weak injection as the signal, while the other with vacuum input as 
idler. When the weak signal is amplified by the strong pump field via the FWM process, 
an idler beam is generated, and the frequency overlap between the amplified signal and 
generated idler fields is negligible for the non-degenerate FOPA. 

When the FOPA is pumped with a single frequency CW laser, as illustrated in Fig. 1(b), 
the two pump photons have the same frequency mode {upi = ujp2 = ujp), a pair of signal 
and idler photons in the frequency modes Us and Ui are uniquely defined. In this case, 
the input-output relation of the single mode FOPA is given by the well-known Bogoliubov 
transformation: 

bs,i = fias,i + i^alsy (1) 

where a and b refer the input and output modes, the subscript s, i respectively denote the 
signal and idler fields, the coefficients /i and v are determined by the parametric gain of 
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FIG. 1: (a) A typical diagram of the non-degenerate phase-insensitive FOPA. Fg and Fj, filters; 
Dl and D2, detectors; Ag and j4j, radio frequency amplifiers; DDS, differential detection system, 
(b) The sketch map of the frequency modes for the single frequency laser pumped FOPA. Two 
pump photons of a single frequency laser uniquely defines the single frequency modes of a pair of 
signal and idler photons, (c) The sketch map of the frequency modes for the pulse pumped FOPA. 
A single frequency signal mode (dotted arrow) will be correlated to many idler frequency modes 
(enveloping curve in gray), and vice versa. Since the pump is in multi-frequency modes, the pulsed 
signal and idler beams are also in multi- frequency modes (enveloping curves in black). 



FWM and satisfy the relation |/ip — I'^P = 1- For input signal in a coherent state \q)s with 
photon number much greater than 1, i.e., 



lar > 1, 



(2) 



it is straightforward to deduce the expressions of all the parameters of the FOPA by calcu- 
lating the quantum average of the corresponding operators over the state \^)'^^^ = |a)s |0)i, 
where |0)j denotes the vacuum input of the idler field and the superscript (s) denotes the 
single frequency mode case. 

The amplified signal (generated idler) beam is then measured by the detector Dl(2) 
after passing through the filter Fs(j^ to reject the strong pump (Fig. 1(a)). The average 
photon number /^(j) and the photon number fluctuation A/^.^^ of the signal (idler) fields are 
converted to the DC and AC currents of Dl and D2, respectively. For identical Dl and 
D2, when the efficiencies of the detector Dl(2) and filter Fs(j) are ideal, the average photon 
numbers of the amplified signal and generated idler beams can be written as 



/W = (6tS^) = |;^|>|2+|z,| 



M - ih^ 



bi) = \i^\^\a\^ + \iy 



(3a) 
(3b) 



where h*^ (i)bs,{i) is the photon number operator of signal (idler) beam. In our theoretical 
model, we assume the condition |ap ^ 1 is always satisfied, so the contribution of sponta- 
neous terms |z/p in the right hand side of Eqs. (I3^) and (jHb) are negligible. According to 
the definition of the photon number gain — the ratio between the average photon numbers 
of amplified signal and input signal 

9 = f , (4) 

-'in 

we have g^'^^ = |/ip. Moreover, according to the definition of the intensity noise of the 
amplified signal and generated idler fields, we have 

A/i^)' = {{b%r) - {b%f = (|/.r + l/il Vnia^ (5a) 

A/f ' = iibM) - mr = (I/.I vp + wna\\ (5b) 

After using the SNL of signal (idler) beam AJ^^^.^^ = (&Lj)^s(j)) to normalize the intensity 
noise AI ,1 , we obtain 

i?(^) = i?f ) = 2^(^) - 1, (6) 

indicating the normalized intensity noise of the amplified signal beam nf is equal to that 
of the generated idler beam R^ . Notice that they are both above the SNL because of 
amplification. 



As for the noise figure (NF) of the FOPA, according to the definition 10 | 



with 



m 

and 

I' 



SNRout 



SNR,n = ^ = l«l' (8) 



SNRout = -TTJ (9) 

A/i 

respectively denoting the signal to noise ratio of the input and amplified signal beams, the 
NF of the single mode FOPA is given by 



\i^ 



2 



2^(*) - 1 



iVF(«) = 1 + f^ = -^^^, (10) 

\H 9 

which shows that the parametric amplification process adds the excess noise to the amplified 

signal beam and the 3 dB quantum limit of NF can be achieved in the high gain limit. 



In the ideal case, the photocurrents of Dl and D2 are highly correlated and can be 

A(f _f-)'2 

subtracted from each other to yield quantum- noise reduction, i.e., . .Vi Af^ < 1. In practice, 
the measured noise reduction is very sensitive to the collection efficiency of twin beams 
and the response function of differential detection system (DDS), consisting of Dl and D2 
followed by the radio frequency (RF) amplifiers Ag and Ai, and the subtractor (see Fig. 
1(a)). For brevity, we assume throughout this paper that the transmission efficiencies of the 
filters Fg and Fi at their central wavelengths are perfect. In this situation, the collection 
efficiencies of the signal and idler beams at the single frequency Us and Ui are perfect, so we 
only need to analyze the influence of DDS with non-ideal response. 

The response function of the DDS is described by Qi{^) and Q2{^), which are determined 
by the quantum efficiencies and electrical circuit of Dl and D2, respectively. Here Q denotes 
the radio frequency. The DC response of the DDS, Qi(2){0), is associated with the quantum 
efficiency through the relation (5i(2)(0) = ^^'^ where e is the electron charge, r]s(i) is the 
quantum efficiency of Dl (D2) , and h is the reduced Planck constant; while the AC response, 
<5i(2)(^) (^ 7^ 0), is not only determined by the response of Dl (D2) Di(2)(f2), but also 
determined by that of the RF amplifier As(^i){Q). For convenience, the ratio r between the 
AC response Qi{Q) and Q2{^), defined as 

which can be adjusted by changing the electronic gain of the amplifier As or Ai. 

In practice, the quantum efficiencies of Dl and D2 are not perfect, i.e., f]s < I and rji < 1. 
Modeling the non-ideal detectors as the beam splitters, which couple the vacuum mode Vg 
and Vi to the signal and idler fields, the field operator incident on Dl and D2 are given by 



Cs = y/risbs + a/I - VsVs, (12a) 

Ci = ^/rj'ibi + a/1 - TjiVi. (12b) 

In this situation, the measured photon numbers of the amplified signal and idler beams in 

Eq. (3) are rewritten as 

lf = {c%)=r^sWW? + \^?) (13a) 

/f = (clQ)=r^.(|H>P+IH'), (13b) 

and the intensity noise of individual signal and idler beams in Eq. (5) is rewritten as 

Hlf? = {{c%f) - {{c%)Y = {2r^lHM' + Vs\^?)W\\ (14a) 



A(/i^^? = {{c\c,f) - ((cIq))2 = {2v^\u\' + r^M')H'- (14b) 

For the operator of the measured intensity difference of the twin beams 

f,'^' = c% - r4c., (15) 

the normahzed intensity difference noise of the twin beams is expressed as 

A(/r)2 + r2A(/(f)2 A(/(f)2 + r2A(/(f)2 ' 

where A(/^^L)^ = Ijl is the corresponding SNL of the detected signal (idler) beam and 
Hsi = 2?7s?7j|/ip|z/p|Q;p is the quantum correlation term. Substituting Eqs. (IT3l)-(IT^ into 
Eq. f lT^ . we obtain 



I 1 9 I I 1 9 9 



Rt {Vs,Vi,r) = ^ ^ . (17) 



To illustrate the factors influencing the value of B!^ {r]s,r]i,r), we first plot -R^ as a 
function of the gain g^^^ by assuming the ratio r in Eq. (11) is equal to 1 and by varying 
the efficiencies rjs and r^j, as shown in Fig. 2(a). When the efficiencies of both Dl and D2 
are 75% or 85%, R^ decreases with the increase of g'^^'^] for a given g^^\ R^ obtained for 
Vs = Vi = 85% is lower than that for i]s = rji = 75%. When efficiencies of Dl and D2 are 
rjs = 75% and rji = 85%, respectively, R[ is higher than that for 7]s = 7]i = 75% in the 
high gain regime, however, R^ can be slightly lower than that for r]s = r]i = 85% when the 
gain is within the regime of 3 < g^^^ < 10. Therefore, for the case of rjg = rji, R^ decreases 
with g^^^ and increases with the decrease of i]s(i); while for the case of r]s 7^ i]i, R^ does 
not always decrease with g^'^^: after obtaining the minimum R[ at a certain gain, R^ will 
start to increase with g^'^\ 

Eq. f lT7|) shows that for the given values of t]s, rji and g^^\ R[ can be minimized when 
the ratio r takes the optimized value 






We then plot R^ as a function of g^^'^ by varying the efficiencies 77^ and rji for the case of 
r = r^pj. As shown in Fig. 2(b), no matter rjs and rji are equal or not, R[ always decreases 
with the increase of g^^\ Moreover, for a fixed value of g^'^\ our calculation indicates that 
R[ decreases with the increase of r]s and r/j. 
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FIG. 2: (a) The noise reduction R^ as a function of the gain g^^' under different quantum 
efficiencies for the cases of (a) r = 1 and (b) r = r^J^. 

We note that it is impossible to obtain the ideal noise reduction R^ {r]s,r]i,r) = for 
the input signal with photon number |ap 7^ 0, even if the detectors Dl and D2 are perfect. 
To clearly illustrate this point, let's look at the expressions of R^ when the condition 



Vs 



rji = 1 is fulfilled. For the case of r 



R 



l{s), 



i,i;'^opt) 



'"opj, Eq. flT7|) is rewritten as 
1 



(H + IH)" 

while for the case of r = 1 , Eq. ( iT7|) is rewritten as 



(19) 
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(20) 
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l/^l + r I 

Both Eqs. flTOj) and fl20l) show the intensity difference noise of twin beams is always less 
than the SNL (1 or dB), however, the condition R^ > is always satisfied. Moreover, 
a comparison between Eqs. f ll9|) and fl2U]) indicates that in the high gain limit, the noise 
reduction can be improved up to 3 dB by optimizing the ratio r. 



III. QUANTUM DESCRIPTION OF A PULSE PUMPED FOPA 



When the FOPA is pumped by a strong pulsed field having multiple frequency com- 
ponents, such as a mode-locked laser with a pulse width of a few pico-seconds, a single 
frequency signal mode is correlated to many idler frequency modes, and vice versa (see the 
dotted arrow and the corresponding enveloping curve in Fig.l (c)). In this situation, both 
amplified signal and generated idler beams are pulsed fields with broad bandwidth, as illus- 
trated by the spectrally broadened curves in Fig. 1(c). Therefore, it is necessary to deal 



with the multi-frequency modes. In this section, we first extend the single mode Bogoli- 
ubov transformation to the multi-frequency domain. Then, according to the Hamiltonian of 
FWM in fiber, we deduce the four transformation functions of the multi-mode Bogoliubov 
transformation, from which the general expressions of the key parameters of a pulse pumped 
high gain FOPA can be derived. 

A. Multi-mode model of a pulse pumped FOPA 

An optical field propagating along the fiber (denoted as the z direction) can be quantized 
by using one dimensional approximation if the polarization mode is well defined. Therefore, 
the pulsed field operator in optical fiber is written as [25|] 

Eit)=E^+\t) + E^-\t), (21) 

where the positive frequency field operator £■(+) (t) and negative frequency operator E^^'' (t) 
satisfy the relation E'^+^t) = [E'^-\t)]\ and E^+\t) is given by 

^(+)(t) = ^ f dua{u)e'^^-'-'''^ = -^ [ dua{u)e-'^''. (22) 

V 27r J V 27r J 



Here the annihilation operator at frequency u, a{u), satisfies the commutation relation 
[a{co) , a"^ {uj')] = 6{u — u'). Note that t' = t — kz/u, so the distance in propagation is 
equivalent to a delay in time. 

For the pulse pumped FOPA, the strong pump pulses remain classical, but the signal 
and idler fields are quantized. Extending Bogoliubov transformation in Eq. (jT) to the 



multi-frequency domain and using Eq. (22), we write the input-output relation as 23| : 



bs{ujs) = U^ctsiuj'sW = / his{uJs,uj'^)ds{uj'^)du'^+ I h2s{uJs,^'i)a\{uj'i)duj[ (23a) 

hi{uji) = U^ di{uj'i)U = j hu{uji,uj[)di{uj[)duj[+ / h2i{uji,uj'^)d\{uj'^)duj'^, (23b) 

where the unitary evolution operator U is determined by the Hamiltonian of the FOPA (see 
later for the explicit form), the footnotes S*, / denote the frequency ranges of the signal and 
idler fields, and the four Green functions hnj{co,u') (n = 1,2 j = s,i) are referred to as the 
transformation functions. To ensure the satisfaction of the commutation relations of the 
operators ds(i){u'^r^\) and &s(i)(ws(j)), the transformation functions should be constrained by 

10 



the relations: 

hu{us,uj'^)h2i{ui,u'^)duj'^- h2siuJs,(^'i)hii{ui,u-)dUi = (24a) 

Ji 

his{uJsl,Uj's)hl,{u}s2, Uj's)duj's - / h2s{u} si, Uj'i)h2,{uJ s2, ^'i)duj'i = 6{uJsl - UJs2) (24b) 



hii{uii, u'i)hli{uji2, uj'i)du'i - / h2i{ujii, uj'^)h*2i{uji2, uj's)duj'^ = 6{uii - Ui2). (24c) 
/ Js 

The weak input signal pulses, synchronized with the pump pulses, is ideally in a multi- 
mode coherent state 

\a)s = \{a{ujs)}) =exp\ / [a ■ s{uJs)d'^ - h.c.]dus\\0), (25) 



where a is the complex amplitude, and the frequency distribution function s{cOs) satisfies the 
normalization condition J \s{us)\'^dus = 1 so that the photon number of the input signal 
is lin = |ap, which is the same as Eq. ([2]). In general, the bandwidth of the input signal is 
much smaller than that of the central frequency. So the quasi-monochromatic approximation 
applies, and the frequency integral range from to oo can be treated as from — oo to cxd. 
For the sake of brevity, all the integral ranges from — oo to oo will be omitted hereinafter. 

In contrast to the single mode FOPA, wherein the filters Fg and Fi with ideal transmission 
efficiency at the central frequencies result in perfect collection efficiency of twin beams, the 
collection efficiency of the broad band multi-mode twin beams is associated with the spectra 
of Fs and Fi. In this case, the field operators incident on the detectors Dl and D2 are given 
by 

Csiujs) = ^/v'sfs{(^s)bsiuJs) + ^a/I - r]J^{us)vsi(^s) (26a) 



Ci{uJi) = y/rhfi{^i)Hu}i) + i\j^ - riif'l{ui)vi{uji), (26b) 

where the complex function fj{ujj) {j = s,i) describes the spectrum of the filter Fj. Eq. 
( l26l) implies that the measured quantum noise property of the pulse pumped FOPA is not 
only influenced by the quantum efficiencies of Dl and D2, but also influenced by the spectral 
functions fs{oJs) and fiiuji). For the ideal detectors, we should have r/^ = r/, = 1, while for 
the perfect collection of the twin beams, we should have \fs{uJs)\ = \fi{uJi)\ = 1. 

Similar to the single mode FOPA, the key parameters of the pulse pumped FOPA can be 
obtained by calculating the quantum average of the corresponding operators over the state 
|\E') = |a)s|0)i. In the next subsection, we will derive the analytical forms of hnj{ou,u') {n = 

11 



1,2 j = s, i), which are the key for deriving the formulas of the operators of twin beams (see 
Eq. (ESD). 



B. Derivation of the transformation function 



The unitary operator in Eq. (!23|) is determined by the Hamiltonian of parametric process. 



So let's begin with t 



optical fibers [26 



le Hamiltonian of the co-polarized pulse pumped FWM in single-mode 



27| 



Hit) = C^x^'^ J dV[E,,it)E,2it)Ei-\t)El-\t) + h.c], (27) 

where Ci is a constant determined by experimental details and the units of quantized optical 
fields, x^^^ is the 3rd-order real nonlinearity, Epi and Ep2 are the classical pump fields, 

Ei-\t) = -^ j duj'M^',)e-<'^^-<'^ (28a) 



and 

E\-\t) = -^ f duj[a\{uj[)e-^'~'^^^-<'^ (28b) 

y 2,71 J 

are the quantized negative-frequency field operator of the signal and idler beams, respec- 
tively. The Gaussian shaped strong pump pulses propagating along the fiber can be written 



as 



27 



Epn{t) = Eoe-'^^^' / e-Kn-^po)V2-|g»{fc^„.-a;p„i)^^^^ ^^ _ ^^2), (29) 



where Up, Upo and kp are the bandwidth, central frequency, and wave vector of the pump, 
respectively, and Eq is related to the peak power through the relation Pp = 2tt(TpEq. The 
nonlinear coefficient 7 is expressed as 7 = -^^ — , where Aejf denotes the effective mode 
area and c is the speed of light in vacuum. After substituting Eqs. ( 128|) and fl29|) into Eq. 
f l27|) and changing dV in Eq. fl27|) to dV = A^fjdz, we carry out the integral over the whole 
fiber length L (from -L/2 to L/2) and arrive at the Hamiltonian in the time dependent form 

^f^*^) = — I — '' 2 2^^ I dujpidujp2duj'^duj[dl{uj'Xa\{uj[) 

2 ^ 2(7^ 

where A/c = kg + ki — 2kp + 27Pp is the phase mismatching term, and the term 27Pp is 
originated from the self-phase modulation of pump. 

12 



To obtain the formula of the unitary evolution operator U = exp {- — jj-^}, we need to 
carry out the integration over the time and over all the possible combinations of Upi and 
Up2 within the pump bandwidth. Since the time integral gives rise to the 6 function in 
Auj = Upi + Up2 — oj'g — uj[ to guarantee the energy conservation at single-photon level, we 
arrive at 

f/ = exp {g[ I U{J^,uj',)al{uj',)a\{uj[)duj'M - h.c]], (31) 



where the so called two-photon joint spectral function (JSF) of FWM has the form of 291] 



^K,c.J = ^^7^exp{ — }.^nc(^), (32) 

which is determined by the pump envelop function exp { 4^2 '^^° } and phase matching 

function sinc{^^), and is referred to as the probability amplitude of simultaneously finding 
a pair of signal and idler photons within the frequency range of wl, — )■ w^ + dco'^ and u'^ — )■ 

—SiCi'yP LcA^ 

uj[ + dijj[, respectively. The coefficient G = 3^^^ ^ ^^^ in Eq. f l3T]) determines the gain 

of FWM, and the coefficient C in Eq. ( 132|) is a constant used to ensure the satisfaction of 
the normalization condition // \ip{Ug,Lo'j)\'^dLOgdco'^ = 1. According to Eq. ( !32|) . the JSF is 
generally asymmetry since Ak is frequency dependent. Therefore, 'ip{ui,U2) = 1^(^^02,001) is 
generally not satisfied. 

To obtain the formula of the transformation functions, we rewrite the input-output rela- 
tion by substituting Eqs. ( 131]) into Eqs.( l23|) and using Baker-Hausdorff lemma: 

bjicoj) = Wdj{co'j)U = dj{u'^) + [djioo'^), B] 

+ ^[«.K-)> [«.K-),^]] + ^[«.K), [«.K-), [aA^-),m + ... , (33) 

where the operator B is expressed as 

B = G [[[i,{u',,u[)d\{u',)d\{u[)duj'M - h.c]. (34) 



Then the transformation functions can be written in the form of an infinite series 
his{uJs,uj',) = 6{uJs -uj'^) + ^ ■—- / / ■ • ■ / duJidu2 ■ ■ ■ duJ2n-l 

[iIj* {u2,uji)ij* [ooA^i^^iW {oo&,(^b) ■ ■ ■ i'*{(^'s,(^2n^i)\} (35a) 
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h2s{uJs, W-) = Gip{ujs, W-) + ^ . ., / / • • • / duiduj2 ■ ■ ■ duJ2n 

{ [i'*{uj2, UJi)lp*{uJi, W3) ■ ■ ■ V'*(w2n, UJ2n-l)] 

[ipiUs, Ui)lp{uj2, UJs)lp{uj4, W5) ■ ■ ■ V'(w2n, W-)]} (35b) 

hu{uJi, Uj'i) = 6{uJi - Lo'i) + ^ — — y //■■■/ C?WlC?W2 " " " duJ2n-l 

{[llj{Ui,Ui)lp{uj3,U2)ij{0J5,U4) ■ ■ ■ 1p{uj2n-l, UJ2n^2)] 
[lp*{uJi,UJ2)'lp*{uJ3,U4,)lp*{uj5,UjQ) ■ ■ ■ V'*(w2n-1, W-)] } (35c) 

h2i{uji, u',) = Gij{uj',, Ui) + ^ / / ■ ■ ■ / dujiduj2 ■ ■ ■ duJ2n 

{ [i)*{Ui, UJ2)lp*{0J3, U4) ■ ■ ■ 1p*{uJ2n-l,UJ2n)] 

[lp{uJi, U-i^lpiUs, ^2)^(^5, ^4) ■ ■ ■ i'iuj's, W2n)]}- (35d) 

Eq. ( l35|l clearly shows that the transformation functions satisfy the relations 

his{uJs, uj's) = h\^{u[, Us) (36a) 

hii{uuu[) = hli{Ui,Ui) (36b) 

h2siuJs,(^'j) = h2i{uJi,UJs), (36c) 

indicating that there is a correlation between the signal and idler beams. 

It is worth noting that Bogoliubov transformation in multi-frequency modes (Eq fl23l) ) 
with the transformation function shown in Eq. (!35|) is also suitable for describing the single 
mode FOPA. For example, if the JSF in Eq. fl32|) takes the limit of a single frequency laser 
pumped case, i.e., ip{Ug,u'j) = 5((u;^ ~ ^so) + {^i — ^io)), with Uso and Wjo denoting the 
single frequencies of the signal and idler beams, the input-output relation of Eq ( |23|) would 
transform into bs(i) = ^ds(i) + J^d},. with yU = coshG and u = sinhG, which is exactly the 
same as Eq. ([1]). Hence, Eqs. (!23|) and (!35|) can be viewed as a general description of the 
input-output relation of an FOPA. 

C. Key parameters of the pulse pumped FOPA 

For the measurements of the signal and idler fields, in addition to taking the quantum 
efficiencies, tjs and r/j, and the spectrum of filter ^^(j) into account, we need to consider the 
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response times of the detectors Dl and D2 as well. When the response time of each detector 
is much longer than the pump pulse duration, but is much shorter than the period of two 
adjacent pump pulses, the average photon number of the amplified signal (generated idler) 
beam per pulse can be expressed as 



Isii) = {Isi^)) = I {Ws(^){m)dt, (37) 

where 



^(-)m7?(+) 



I,{t) = Ep{t)Ep{t) (38) 



with 



if )(t) = -= / c,ico,)e<'^^--^'^duj (j = s,t) (39) 

v2vr J 

is the photon number operator of the pulsed field. The time integral ranges in Eqs. ( 137|) 
and ( l39l) are omitted, because they can be treated as from — oo to oo. 

With the filter operator incident on Dl and D2 (Eq. ([26])), we substitute Eqs. ([23]), fl25ll . 
( I39l) into Eq. ( 1371) . and obtain the expressions of the detected photon number 



= |a| r]s / dujs\fsi(^s)\ SsiuJs) 

+ 7]s dus duj'^\fs{us)\^\h2s{uJs,uj'i)\^ (40a) 



= |a| r]i / dui\fi{ui)\ Si{ui) 

+ rii f du, [ du',\fi{u,)\^\h2M.^'s)\\ (40b) 



where 



and 



Ss{^s)= hl^{ujs,uj'si)hu{ujs,uj's2)s*{(^'si)s{uj',2)duj'^iduj'^2 (41a) 



are the power spectra of the signal and idler beams. The first terms in the right hand side of 
Eqs. (l40k ) and (140b) are originated from the stimulated emission of input signal, while the 
second terms are from the spontaneous emission. Since we have assumed that the photon 
number of the weak input signal per pulse is much greater than 1, the spontaneous emission 
terms are negligible. 
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Using Eq. (140 aP and the definition of plioton number gain in Eq. (jlj), we find the 
expression of the photon number gain of the pulse pumped FOPA is 

g = f Ss{cOs)dLo, (42) 

for the case of r^^ = 1 and \fs{ujs)\ = 1- 

To compute the intensity noise of the amphfied signal (generated idler) beam 



A/,^(,) = (/.(.)/.(.)) - /i(,) = jj {Wsi^){tl)Is(^){t2)\^)dt^dt2 - /i(,), (43) 

we substitute Eqs. ([23]), ([25]), ([26]) and ([37])- ([39]) into Eq. ([43]), and arrive at the detailed 
expression 

A/2 = {H,, + H,2 + H,s) |«P (44a) 

Alf = {Ha + Hi2 + H,i) \a\\ (44b) 



with 



Hsi= vl 1 1 ■■ ■ j dusidus2duj'gidu'^2d'(^s3 

{\fs{uJsl)\'^\fs{uJs2WK^{uJsl,Uj'si)his{uJs2,(^'sl) 

Ksi^sl, Uj'^2)s*i^'s2)hls{^s2, ^^3)5(^^3)} (45a) 

Hs2= ^^s ■■■ dujsidujs2duj-^duj'^^duj'^2 

{\fs{uJsl)\'^\fs{uJs2)\'^K^{uJsl,UJa)h2s{(^s2,(^il) 

Hii= r]f ... dujiiduji2duj'^^duj'^2(^uj'^3 

{\fi{(^il)\'^\fi{uJi2)\'^h*2i{uJiuUj'^^)h2i{uJi2,Uj'^l) 

h2i{ujii,uj',2)s*{uj',2)h*2i{uji2,u',^)s{u'^r^)} (45c) 

Hi2 = ^^ ■ ■ ■ duiidui2du[idu'g^du'g2 

{\fi{.^il)?\fi{.^i2)?hii{uJii,Uj[^)hl^{uJi2,Uj[^) 

h2i{uJi2, uj',,)f*{uj',,)hUujiuUj',2)f{uj's2)} (45d) 
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and 



Hvs= Vs I dusil-r]s\fsi(^s)\ ) 

d^'siM2[Ksi^s, w^i)/ii.(w., uj',2)s*{uj'si)s{uj',2)] (46a) 



Accordingly, the normalized intensity noise of signal (idler) beam is 

Moreover, we have SNRout = # = \'^\'''i^sUdu^^\fs(u^^?Ss{us)? ^ substituting Eqs. gOl) and 
P7|) into Eq. (|9]). According to the definition in Eq. ([7]), we obtain the formula of noise 

figure 

NF = ^^1 + Hs2 + H,s . . 

^lUdus\fs{usWS,{uj,)Y- ^ ^ 

Eqs. (H^ - (H7|) indicate that the intensity noises of the individual beams AJ| are not only 
related to detection efficiencies, rjs and 77,, but also to the spectral functions of filters, fs{oJs) 
and fiiui). When 77^ and rji are very low, i.e., rjj ^ 1 (j = 5,2), the intensity fiuctuation 
^llu) in Eq. ( l44l) is dominated by the term if^j (j = s, i), while the contribution of the term 
Hjk {k = 1,2 j = s, i) is negligible. On the other hand, when the bandwidth of the filters 
Fg and Fi are much narrower than that of input signal, A/^/^^ is also dominated by the term 
Hyj {j = s,i). Since the non-ideal detector induces the vacuum noise to the individual signal 
or idler field, the former is easy to understand. However, the understanding of the latter is 
not so straightforward. We think this is because, when the signal and idler beams are in 
multi-frequency mode, the filter Fg {Fi) with bandwidth narrower than that of the amplified 
signal (generated idler) beam can also be viewed as a loss, which introduce vacuum noise to 
the signal (idler) beam as well. Hence, the term H^j (j = s, i) in Eqs. fl46ap and fl46bp is 



originated from the loss induced vacuum. If the detection and collection efficiencies of the 
signal and idler beam are ideal, we would have H^j = (j = s,i). 

Finally, we derive the expression of the intensity difference noise of the twin beams by 
using the definition 



Mt = {It) - {hY = 1 1 {Wt{ti)it{t2mdt,dh - { I {Wt{m)dty, (49) 
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and 



where the operator It = Ig — rli is the photon number difference with a weight factor of AC 
response ratio r. Using Eqs. (|23|) . ( 126|) . fl37|) -f l39|) and P3l) . Eq. (H9|) is transformed into 

AI^ = [{Us) + r^{hh) - 2r{iM - (J, - rh)'' 

= AP^+r'Alf-2r{H,a + H,,2)\a\\ (50) 

where 

Hsii = VsVi / / • • • / dusduiduj'g^duj'g^d^'ss 

Ksi^s, u's2)s*{u's2)hli{uji, uj's3)s{u's3)} (51) 

Hsi2 = VsVi / / • • • / dusduiduldu'sidu's2 

hisius, u'si)f{uj'si)h2i{uji, uj's2)r{^'s2)} (52) 

are positive terms originated from the quantum correlation between signal and idler twin 
beams. It is obvious that Alf is less than the sum of intensity noise of individual beams 
{Igls) +f'^{lili) due to the correlation of twin beams. Consequently, similar to Eq. ( TT6|) . the 
general expression of the normalized intensity difference noise for pulsed twin beams 

_ [Hsi + Hs2 + H,s + r^{Hii + Hi2 + H,i) - 2r{Hsii + Hsi2)]\a\^ ._„s 

is obtained by substituting Eqs. (jH]) and fl5T]) - fl52|) into Eq. ( 150|) . In addition, Rt can be 
further minimized if the AC response ratio r takes the optimized value 

^4Ui{Hsii + Hs^2? + [Is{Hii + Ha + H,i) - h{Hsi + Hs2 + H,s)f 



f — fopt 



'^hi.Hsii + Hsi2) 
Ii{Hsi + Hs2 + H,„) - Is{Ha + Ha + Hvi) 



(54) 



2Ii{Hsii + Hsi2) 

From the formulas of the photon number gain, noise figure and intensity difference noise 
in Eqs. ( H2l) . ( H8|) and ( l53l) . one sees that the key parameters of the FOPA are determined 
by the transformation functions in Eq. ( I35l) . which highly depend on the JSF in Eq. (1321) . 
In the following sections, we will focus on studying the noise characteristics of the pulse 
pumped FOPAs with two typical kinds of JSF. 
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IV. PULSE PUMPED FOPA WITH SPECTRALLY FACTORABLE JSF: 



We first study the FOPA with a spectrally factorable JSF. Under this condition, both the 
signal and idler modes are in single temporal mode, so we expect to recover the results for 
single mode FOPA in Sec. |IT1 For the sake of brevity, in this section, we assume the quantum 
efficiencies of Dl and D2 are perfect, and the collection efficiency of the twin beams is ideal. 

A factorable JSF, obtained by properly regulating the pump pulses and tailoring the 



dispersion of optical fibers, is written as 29 1 



i;^S\uJ'^M) = <t>s{0*V^H), (55) 

where 0s(w!,) and V'j(wQ with the normalization condition of j \(f)s{uj'g)\^duj'g = 1 and 
f \ipi{cOi)\'^dUi = 1 are the probability amplitude of finding a pairs of signal and idler photons 
within the frequency range of w^ — )■ u'^ + dUg and uj[ — )■ u^ + du'^, respectively. In fact, looking 
at Eq. (!55l) from another point of view, we find (psi^'s) and fi{co'i) can also be viewed as the 
gain spectra in signal and idler fields, respectively. 

With the JSF in Eq. (!55|) . the transformation functions in Eq. (135|) are simplified to 

h[i\uJs,co',) = 5ius - cu'J + (coshG - 1) • (P^ico^WM) (56a) 

hi^Jius^uj'i) = sinhG ■ (f)s{uJs)^i{(^'i) (56b) 

h\i{uji, u;-) = 6{uJi - w-) + (coshG - 1) ■ ipi{uji)ip*{uj'i} (56c) 

h^i {u)i,u'^) = sinhG ■ (pi{ui)(psiu)'s)- (56d) 

Note that in this section, we will use the superscript (/) to represent the result of an FOPA 
with the factorable JSF. 

With the JSF in Eq. ( l55l) . the average photon numbers of the amplified signal and 
generated idler beams per pulse in Eq.( l40l) are simplified to 



jf/) = (l + |F-sinhG'p)|a|2 (57a) 

4^) = (|F■sinhG|2)|a|^ (57b) 

where the coefficient F = J s{uJs)(j)*{u)s)duJs is determined by the matching between the 
spectra of signal injection s{ijj) (see Eq. ( l25l) ) and that of the gain in signal band (ps{co). 
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Therefore, we obtain the photon number gain of the FOPA 

(,(/) = l + |F.sinhG'|2 (58) 

by substituting Eqs. (jST]) and (^ into Eq. (g]). 

To show the factors influencing the noise performance of the FOPA with factorable JSF, 
we first substitute the simplified transformation function in Eq. (13^ into Eqs. (145 1) - (H^ and 
( 15T|) -( !52|) . For the case of rjs = rji = |/s(w)| = |/j(i^)| = 1, the expressions of the terms in 
Eqs. (|15])-(116D and ([5l])-([52]), which determine Alf (j = s,i) in Eq. (HH) and A/^ in Eq. 
(!50|l . are simplified to 

H^P = l + \F\^-{\coshG\^-l), (59) 

hIP = \F\^-\smhG\\ (60) 

Hi2 = h\P = H^il = His = \F- sinh G cosh G|^ (61) 

ifif = 0. (62) 

Consequently, the normalized intensity noise of individual signal and idler beams are 

(^) _ 1 + |F|2(| sinh 6-1^ + I sinh G cosh Gp-l) 
^' ~ l + |F-sinhG|2 ^^^^ 

and 

r\-^^ = I sinhGI^ + I coshG|^ (64) 

respectively. Moreover, according to Eq. ( 157|) . the noise figure in Eq. (HHj) and the nor- 
malized intensity difference noise of twin beams in Eq. ( 153|1 can be respectively simplified 
to 

iVF(/) = ^ (65) 

and 

^"'= 2|F..inhGp+l '-^- («) 

Equations ( I63l) - (l66ll show that except for the normalized intensity noise of idler beam RI , 

the remaining three parameters, W\ NF^^ and Ri'\ are associated with the matching 

coefficient F. For the case of F = 1, if we respectively replace the terms 1 + | sinhGp and 
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IsinhGp with |/ip and |z/p, the analytical expressions of the key parameters (Eqs. (163|1 - 
fl66|) ) will be the same as those of the single frequency mode FOPA in Sec. [TTl The results 
indicate that our multi-mode theory of the pulse pumped FOPA are valid. 

An marked difference between the FOPAs described by single temporal mode and single 
frequency mode is the spectra of signal and idler beams. The spectra of the amplified signal 
and generated idler beams obtained by substituting Eq. (156|) into Eq. (l4T!) are written as 

Si^\us) = \s{us)\^ + \F ■ 0.K)n coshG - 1|' 

+F*(cosh G - l)(PsiuJs)(p:icOs) + F(cosh G - iy^:{cUs)M^s) (67a) 

and 

S^/\co,)= \smhG\^\v,{u,)\^ (67b) 

respectively. From Eq. (1^71) . one sees that the spectrum of generated idler is always the 
same as the spectrum of the gain in idler field, however, different from the FOPA pumped 
with a single frequency laser, the spectra of the injected signal and amplified signal beams 
are different unless the condition F = 1 is fulfilled. 

To further understand the difference between the FOPAs described by single temporal 
mode and single frequency mode, we then compare the noise performance of the two cases. 
We find that for a given value of G (or g), the noise figure and the intensity difference noise 
of the twin beams in Eqs. f lB5]) and fl^^ are worse than that predicted by Eqs fITU]) and fl2D]) 
unless F = 1. In practice, it is very difficult to maintain the spectral matching condition 
(psi^s) = s{us) in the high gain regime due to the self- phase modulation and cross- phase 
modulation induced spectral broadening [30(]. Hence, it is very challenging to experimentally 
reahze a pulse pumped FOPA capable of being described by a single mode theory. 

V. PULSE PUMPED FOPA WITH SPECTRALLY NON-FACTORABLE JSF: 



For the conventional optical communication systems, an FOPA with an extremely broad 
gain bandwidth is desirable. When the gain bandwidth of the FOPA is much broader than 
the bandwidths of the pump and input signal, the JSF can be simplified by assuming the 
perfect phase matching condition AA; = 0. In this case, the JSF in Eq. (I32l) is rewritten 
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as [28| 

2vvrcrp 4cr^ 

which is obviously non- factorable. In this section, we focus on numerical investigation of 

the quantum noise performance of this kind of FOPA. 

With the JSF in Eq. ( |68l) . the transformation functions in Eq.( !35l) are reformulated as 

= 6(uj - w') + V -^= e~^^^^^ (69a) 

= V ^e H<^'^+^^ , (69b) 

where G' = CG , and the superscript (6) represents the result for the FOPA with broad 
gain bandwidth. For the sake of convenience, the frequency variables ujs{i) and wVx used in 
previous sections will be replaced with oj and u' hereinafter. 



A. Key parameters of the FOPA in the ideal conditions 

We first study the factors influencing the key parameters of the FOPA by assuming the 
detectors are perfect, and collection efficiency of the twin beams is ideal, i.e., rjs = rji = 1 
and |/s(w)| = |/j(w)| = 1- Without loss of the generality, we assume the spectral function 
of the pulsed input signal can be described by the Gaussian function 

where a denotes the bandwidth. From Eqs. flB^ and fITU]) . we can rewrite the power spectra 
of the amplified signal and generated idler beams in Eq. (14 ip as 



^ oo oo ,Q,,2ni+2n2 _f^ — U )'^ 



ni=0 n2=0 



-(U -UsoY , / 0-2 



"" '^''P ^2(a2 + 4n2a2) W {a^ + An,a^^)ia^ + 4^2(72) ' ^^^^^ 



00 00 






v^ f^ni^o (2^1 + 1)'(2^2 + 1)! " ^2(a2 + 2(2^1 + 1)^^ 



2^ 



^ ^^P ^2(a2 + 2(2n2 + 1)^2) ^V (a^ + 2(2ni + I)a2)(cr2 + 2(2^2 + 1)^^) ' ^ ^ 
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For all the terms of the infinite series in Eqs. (!7Tb) and (17Tb ). whose order corresponds 
to rii 7^ or n2 7^ 0, their bandwidths are larger than that of the input power spectrum 
|s((X's)p. Moreover, since the bandwidth of the term in infinite series increases with the order 
determined by the integers rii and ^2, the bandwidths of both signal and idler twin beams 
increase with G' . 

Using Eq. (!7T|) . the average photon numbers of amplified signal field and generated idler 
beams in Eq. (HOj) can be rewritten as 



00 00 



Jf) = |«|2 y y i^^L ^ (72a) 



(2ni)!(2r22)!V l + 2(ni + n2)p2 

ni=0n2=0 ^ ' ^ ' \ ^ '■■ 



oo oo 
I'' 

ni=0 712=0 



rW _ l^,|2 \^ \^ l*-^! 

(2ni + l)!(2n2 + 1)! V 1 + 2(ni + n2 + \)v' 






where p = <Jp/a describes the ratio of the pump bandwidth to the input signal bandwidth. 
For all the numerically simulated results presented hereinafter, the coefficient G' is within 
the range of < |G'| < 4, and the series in Eqs. (169|) . ( 17T]) and (172|) are calculated up to the 
order of rii = 10 and n2 = 10. To ensure the correctness of this truncation approximation, 

j(b)_j(b) 

we numerically verify the general relation " ^ ' = 1 (/j„ = |a| ), which is the inherent 
nature of signal and idler twin beams originated from the energy conservation of FWM. 



After calculating the series in Eq. ( 1721) to the order of rii = 10 and n2 = 10, we find that 



j(b)_j{b) 

the difference between the calculated result of " , .a' and "1" is less than 10 , showing 



the validity of the approximation. 

Figures [3]^a) and |3]^b) show the normalized power spectra of the amplified signal and 
generated idler beams. In the plots, Eqs. (17Tk) and (17Tb) are calculated by varying the 
value of G' under the condition a = cTp. As a comparison, we also show the power spectrum 
of input signal |s(c(;)p in Fig. El^a). It is obvious that the bandwidths of both signal and idler 
twin beams are broader than that of the input signal and they increase with G", showing the 
spectrum broadening effect. Moreover, in the low gain regime, the spectrum of amplified 
signal is mainly determined by the input signal, so the bandwidth of generated idler beam, 
which is the convolution of the spectra of the input signal and pump, is greater than that 
of the signal beam. However, in the high gain regime, the bandwidths of signal and idler 
twin beams are about equal because both the spectra of signal and idler beams are the 
convolution of the spectra of pump and its counterparts. 

Figure m plots the dependence of the photon number gain g = /I /|ap upon the pump 
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FIG. 3: The normalized power spectra of (a) amplifed signal and (b) generated idler beams at 
different G' under the condition a = a^ . As a comparison, the power spetrum of input signal 
|s(a;s)P is plotted in (a) as well. 

power for the different ratio p = cTp/cr. In the calculation, G' oc P^ is changed by varying the 
peak pump power P^. One sees that at a fixed peak power Pp, g increases with the decrease 
of p because the overlap of temporal mode, which is required for maximizing the gain of 
FWM, is improved by decreasing the bandwidth of pump. When p is less than 0.1, further 
decreasing the ratio to p — ;■ does not result in an obvious increase in g^ because no more 
space is left for improving the temporal mode overlapping. The result in Fig. H] implies the 
pump with a longer pulse duration (narrower bandwidth) gives a better mode overlapping. 
However, it is worth noting that the pump with a short pulse duration helps to obtain a 
high gain under the condition of low average pump power, which is of practical importance. 
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FIG. 4: Photon number gain g versus the peak pump power for the different ratio "p = Cp/a. In 

the calculation, Ak = 0, r/^ = t/j = \fs{ujs)\ = \fi{uji)\ = 1, the parameter /3 = ^j— — ^^ = 1, 

and G' = fSPp. 



To investigate the quantum noise performance of the FOPA with broad gain bandwidth, 
we first deduce the formulas of terms Hjk (j = s,i, k = 1,2) and H^j (j = s,i), which 
determine AJ| (j = s,i) in Eq. (j44]). By substituting Eqs. ([69]) and ([70]) into Eqs. fliSjl - 
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6]), we obtain 



°° |/^/|2ni+2?i2+2n3+2n4 



^^^^^ (2ni)!(2n2)!(2n3)!(2n4)! 



1 + 2(ni + n2 + n3 + n4)p2 



(73) 



-"s2 - 


- -"i2 




oo 


= 


■ i: 




ni— >n4 



I (^/ 1 2ni+2n2+2n3+2n4+2 



(2ni)!(2n2)!(2n3 + l)!(2n4 + l) 



1 + 2(^1 + ^2 + ^3 + 714 + l)p2 



(74) 



°° |^/|2ni+2n2+2n3+2n4+4 

*' ^ (2ni + l)!(2n2 + l)!(2n3 + l)!(2n4 + l)! 



t (75) 

1 + 2(ni + n2 + ng + n4 + 2)p2 

i7„, = {j = s,t). (76) 

With the above, we then numerically study the dependence of the normalized intensity 
noise of individual beams Rj = AI^/Ij (j = s, i) upon the ratio p. Figure [5](a) plots Rs and 
Ri versus g for the cases of p = 10, p = 1, and p — )■ 0, respectively. From Fig. [5]^a), one sees 
that for a fixed value of g (p), both Rs and Ri increase with p (g). Moreover, for the given 
values of p and g, Ri is greater than Rg. This is different from the results of R,"^ {j = s, i) 
or a single mode FOPA [see Eq. ([6])], but has been experimentally confirmed in references 



9| and |20|. 



We also compute the noise figure NF in Eq. fj48l) . Figure [5](b) plots NF versus g for 
p = 10, p = 1, and p — )■ 0, respectively. One sees that for a fixed ratio p, NF increases 
with g in the low gain regime, and approaches a constant in the high gain regime. While 
for a fixed g, the value of NF increases with the ratio p. It is obvious that for the pulsed 
pump case with p ^ 0, NF in the high gain regime is always greater than the well known 
3dB-limit of a single mode FOPA. 

To study the factors influencing the normalized intensity difference noise of the twin 
beams Rt in Eq. (|53|) . we need to calculate H^^l and H^^l ^^ well. By substituting Eqs. (|69|) 
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FIG. 5: (a) The normalized intensity noise of the individual signal/ idler beams Rj (j = s,i) and 
(b) noise figure NF versus the photon number gain g for the different ratio p = ap/a. In this 
calculation, Ak = 0, and rjs = rn = |/s(f:i^)| = |/i(w)| = 1. 
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FIG. 6: The normalized intensity difference noise of the twin beams Rt versus photon number gain 
g for p = 0.1, p = 1 and p = 10 when r is set to (a) r = 1 and (b) r = ropt, respectively. 



and (170|) into Eqs. (15T!) and (!52|) . we find the following relation 



rib) 



r(b) 






if, 



[b) 



H. 



(b) 



i2 



(77) 



In Figs. [6]^a) and [6](b), we plot Rt as a function of g for different values of p when the AC 
response ratio r takes the values of r = 1 and r = Vopt, respectively. Fig. [6](a) shows Rt is 
independent on p and always decreases with the increase of g; while Fig. El^b) demonstrates 
that Rt still decreases with the increase of g, but Rt can be further decreased and depends 
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FIG. 7: The optimized value rgpt and the photon number ratio Ig /l\ as a function of g. In this 
calculation, the phase mismatching term of the FOPA is A/c = 0, and the detection efficiencies and 
the filters satisfy r]s = rn = \fs{'^)\ = \fi{'^)\ = '^■ 

upon the ratio p. For a fixed gain g, Rt decreases with p, which is different from the case 
of r = 1. Moreover, when p is very small, for example, p = 0.1, the difference between the 
-Rf-values with r = 1 and r = r^pt is about 3 dB in the high gain limit; when p is very big, 
for example, p = 10, i?(-values with r = 1 and r = r^pt are almost the same. 

The results in Figs. |6t^a) and [6](b) not only indicate the noise reduction of Rt can be 
improved by properly reducing the ratio p, but also imply that Vopt might depend upon p 
and Vopt might be very close to 1 when the value of p is greater than 10. This is shown in 
Fig. [TJ in which we plot r^pt as a function of g for p = 0.1, p = 1 and p = 10. In addition, we 
also plot the photon number ratio Is /l\ versus g in Fig. [3 One sees that Vopt is always 
within the range of /i /l\ and 1. For a fixed value of g, Vopt decreases with the increase 
of p. As we expected, Vopt obtained for the case of p = 10 is about 1 if (yf is greater than 3. 
Moreover, in the high gain limit, both Vopt and Is /ll approach 1, which is irrelevant to p. 

B. Influence of the detection loss and collection loss of twin beams on Rt 

In practice, the ideal condition rjs = rji = |/s(w)| = |/j(w)| = 1 can never be fulfilled, 
the measured intensity difference noise Rt is closely related to the loss of the FOPA system. 
So we need to study the influence of the quantum efficiencies {rjs and rji) and collection 
efficiency of twin beams. In this subsection, we assume the filters placed at the output port 
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of the FOPA (see Fig. 1) have a Gaussian shaped spectrum 

/,M = exp{ ~^^^~,^^"^' } iJ = s,^l (78) 

where Ujo and o"/ denote the central frequency and bandwidth of the filter Fj (j = s,i), 
respectively. 

Using the transformation functions, input signal spectrum and the spectrum of Fj {j = 
s, i) in Eqs. fl69|) . flTOj) and fl78|) . respectively, the measured photon numbers of the amplified 
signal and generated idler beams in Eq. (HOl) are rewritten as 

'' ^ ^ {2nM2n2)\ 

ni=0n2=0 ^ ^^ ^ '^' 

{ 5^(4^1^2 + 1) (4^2/ + !) 
+ [2(ni + n2)p' + l] }"'/' (79a) 



^ f2ni + l)!(2n2 + l)! 

ni=0 712=0 

{ s2[(4ni + 2)p2 + l][(4n2 + 2)^^ + 1] 
+ [2(ni + n2 + 1)/ + 1] }-^/^ (79b) 

where s = cr/crf is the bandwidth ratio of the input signal to filter ^^(i). For clarity, we 
add an apostrophe to label the result of a broad-band FOPA with non-ideal collection and 
detection efficiencies. Moreover, by substituting Eqs. fl69|) . flTOj) and flTHj) into Eqs. fj45l) - (jl6j) . 
the expressions of the terms of Alj {j = s,i) in Eq. (14^ are reformulated as 

I /^/ 1 2ni+2n2+2n3+2n4 



00 



i?:f ' = 1', E 



„,^„ (2"i)'(2"2)K2n3)!(2n4)! 



1 + S2ni,2n2,2n3,2n4 



(80) 



°° |/^/|2ni+2n2+2n3+2n4+2 



^^^^J2ni)!(2n2)!(2n3 + l)!(2n4 + l)! 



-I- + S2ni,2n2,2n3+l,2n4+l 



00 



(81) 



I /-" 1 2ni+2n2+2n3+2n4+4 



^^_^^^ (2ni + l)!(2n2 + 1)1(2^3 + l)!(2n4 + 1)! 



^ + S2ni + l,2n2+l,2n3+l,2n4+l 



(82) 
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°° |^/|2ni+2n2+2n3+2n4+l 

^i2 -Vi 2^ 



^^_^^^(2ni + l)!(2n2 + l)!(2n3)!(2n4)! 



l + S2ni + l,2n2+l,2n3,2n4 



^:?^ = ^^EE "" 



(83) 



(2ni)!(2n2)! 



1 I cvac \ 1 _L C^ic 

-■- "T S2ni,2n2,l V S2ni,2n2,2 



°° °° |,^/|2ni+2n2+2 



^ }, (84) 



(2ni + l)!(2n2 + l)! 



y -■- "T S2ni+l,2n2+l,l V '?2ni + l,2n2+l,2 

By substituting Eqs. (EHD, (HOD and ([7HD into Eqs. ([5ID and ([S2]), the other two terms of Rt 
in Eq. (153!) are reformulated as 



|/^n2ni+2n2+2n3+2ra4+2 

(2ni)!(2n2 + l)!(2n3)!(2n4 + l)! 



H'sa = VsVi E 



rii— >ri4 



i (86) 



1 + S2ni,2n2+l,2n3,2n4+l 



°°^ |(-^/|27ii+2n2+2n3+2n4+2 

H'si2 = VsVi 22 



^^^^J2ni + l)!(2n2)!(2n3 + l)!(2n4)! 



-f • (87) 

^ ^ S2ni+l,2n2,2n3+l,2?i4 



Here the coefficients 



^Lmmm = (^1 + k2 + h + h)p'' 

+s^{2 + 4(A;i + k2 + h + h)p^ + Alkiiks + k^) + A;i(2A;2 + /cs + /i;4)]p^} 
+4sy (1 + 2A;i/)(l + 2k2P^){k3 + k^) (88) 

4Tfc2,x = X ■ (fci + k,)p' 

+s^ [1 + 2fcip2] [1 + 2k2v'] (a; = 1, 2) (89) 
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are introduced to simplify the Eqs. (180|) -( 187|) . 

For a fixed value oi p = cTp/a, if the bandwidth of filters -Fs(j) is much narrower than 
that of the input signal, namely, s = cr/o-f — )• C)0, we have ^^ j^^ ^^ ^^ ^ ^fc^fe x- ^^ ^^^^ 
case, the collection efficiency of the twin beams is very low due to the multi-mode nature 



of the signal and idler fields 3l|]- Since the narrow band filters ^^(j) introduce the vacuum 
noise into the individual signal and idler beams, for a given G', the value of H'^j {j = s,i) 
in Eqs fl84|) - fl85|) approaches to I'j (j = s, i) and is dominant over the terms H'-^ {I = 1, 2) 
and H'^^i {I = 1,2) (Eqs. (180D-([83D and ([HID-dHTD). Hence, the normalized intensity noise 
and normalized intensity difference noise approach to the SNL, i.e. Rj — ?■ 1 (j = s, i) and 

Rt-^l. 

In general, the collection efficiency of the twin beams is not only associated with the 
ratio s, but also depends on the ratio p. To reduce the detrimental effect of the collection 
efficiency on Rt, both s and p should be small enough unless the condition s — ;■ or p — ;■ 



is satisfied [3l||. To illustrate this point, we plot the contour of the normalized intensity 
difference noise Rt as a function of p and s in Fig. [H](a). In the calculation, we assume 
G' = 3, r = 1, and the quantum efficiencies of Dl and D2 are perfect. For each combination 
of s and p, we compute the terms in Eqs. fl79|) - fl87j) and substitute the results into Eq. fl53l) . 
Fig. [8](a) shows that Rt depends on both s and p. We notice that with the increase of p, 
the dependence of Rt upon s becomes stronger. When p is less than ~ 1.75, at a fixed p, Rt 
increases with s, but the noise reduction condition of i?( < 1 is always achievable. However, 
when p is greater than ~ 1.75, to obtain the noise reduction of Rt < 1, the value of s should 
be smaller than a certain value, which decreases with the increase oi p. In real experiment, 
it is necessary to use the filter ^^(j) with a certain bandwidth to prevent the strong pump 
background from reaching the detectors Dl and D2. So, to improve the noise reduction of 
twin beams, the pump with a pulse duration longer than that of the input signal, i.e., p < 1, 
is desirable. 

Figure M^a) also shows that the contour Rt = 1 divides the plot into two zones according 
to Rt < 1 OT Rt > 1- For the former case, there is a one-to-one correspondence between the 
combination of s, p, and Rt] while for the latter case, each Rt maps two values of s at a fixed 
value oi p: one is very close to 0, the other is larger. To better understand Fig. [8](a), we 
also plot the contours of the sums of the positive terms and negative terms in Eq. fl53|) in 
Figs. [H](b) and[H!^c), which respectively illustrate the infiuence of collection efficiency upon 
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FIG. 8: Contour plot of (a) the normalized intensity difference noise Rt, (b) (A/^ + Alf) / (Ig + !{)■, 
and (c) 2{Hsii+Hsi2)\oi.\'^ /{Is+h) as a function of the ratios p and s. In this calculation, rjs = rji = 1, 



G' = 3, and r = 1, and Rt 
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FIG. 9: The normalized intensity difference noise of the twin beams Rt versus g for (a) p = 0.1, 
(b) p = 1, respectively. For each setting of p, s in Eg. ([79]) and ([80|) - (j85|) is s = 0.1, s = 1, and 
s = 10, respectively. In this calculation, rjs = r]i = 1. 

the intensity noise of individual beams and upon the correlation of twin beams. For both 
Figs. [8]^b) and[8]^c), the contour decreases with the increase of s (p) at a fixed value of p 
(s). However, since the correlation of twin beams highly relies on the collection efficiency, 
the dependence of the contour on s and p in Fig. |8](c) is stronger than that in Fig. Mjo). 
In particular, when p increases, the value of 2{Hsii + Hsi2)\a\'^/{Is + h) in Fig. Wis) rapidly 
decreases with the increase of s, which is responsible for the multivalued mappings in the 
zone of -Rt > 1 (see Fig. [8](a)). 

To further demonstrate the dependence of the noise reduction of Rt upon the photon 
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FIG. 10: The normalized intensity difference noise of the twin beams Rt versus g for different 
setting of rjs and rji. For the plots in (a), r = 1; while for the plots in (b), r is optimized for 
minimizing the measured Rt- In this calculation, p = 1 and s = 1. 

number gain g, we plot Rt as a function of (7 for s = 0.1, s = 1 and s = 10. In Figs. M^a) 
and[9](b), we assume rjs = rji = 1 and r = 1, and p is fixed at p = 0.1 and p = 1, respectively. 
For each combination of s and p, we compute the terms in Eqs. fl79l) - fl85l) by changing G' 
(g) and substitute the results into Eq. (l53|) . It is obvious that for p and s with fixed values, 
Rt decreases with the increase of g. Moreover, at a fixed value of g and s (p), Rt increases 
with the increase of p (s) due to the decreased collection efficiency of twin beams. 

Having understood the influence of the collection efficiency of twin beams on Rt, we then 
study the influence of the detection efficiency on Rt when the values of p and s are fixed. 
Without loss of generality, we assume p = 1 and s = 1. When the AC response ratio is 
set to r = 1, we plot Rt as a function of g by varying rjg and rji, as shown in Fig. [TO](a). 
It is obvious that for the case of rjs = Vi-, Rt decreases with the increase of rjs^i) at a fixed 
value of g. However, for the case of rjs 7^ rji, Rt is not a monotonically decreasing function 
of g and rjgi^iy Depending on g, the correlation of photon currents originated from the 
quantum correlation of the signal and idler beams might become weak or strong due to the 
unbalanced detection efficiencies. As shown in Fig. fTOT a). Rt for the case of rjs = 75% and 
T]i = 85% is lower than that for rjs = rji = 85% when g is less than 18; however, Rt becomes 
even higher than that for r]s = rji = 75% when g is greater than 30. Making comparison 
between Fig. [TOlfa) and Fig. 2(a), which are obtained for FOPA with non-factorable JSF 
and with single frequency pump, respectively, we find there are two differences: (i) for g, 
rjs and r/j with certain values, Rt in Fig. [TOT a) is larger than that in Fig. 2(a) because the 
collection efficiency of the pulsed twin beams is smaller; (ii) the difference between the curve 
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for Tjs = 75% and rji = 85% and that for rjs = rji = 85% {rjs = rji = 75%) in Fig. [TOlfa) is more 
apparent than that in Fig. 2(a). Therefore, for the pulse pumped FOPA, the measured Rt 



can be minimized by properly adjusting the efficiency ris or r^j [20|]. However, this method 
will introduce extra loss, and the minimized Rt can only be obtained for g with a specified 
value. 

If the ratio r is adjustable, we can minimize Rt by optimizing the value of r without 
introducing the vacuum noise. In this case, as shown in Fig. [TOlfb). for fixed values of r/s 
and rji, Rt always decreases with the increase of g; while for a fixed value of g, Rt always 
decreases with the increase of ris and rji. Hence, for both the signal and idler channels, the 
higher the efficiencies are, the better the noise reduction is. 

C. Influence of excess noise of input signal 

Finally, we analyze the infiuence of the excess noise of input signal on Rt , since the noise 
of input signal pulses originated from the mode-locked fiber lasers is often higher than the 
SNL [20] . When the excess noise is included, the input signal can be viewed as a mixture of 
coherent states, whose density operator and average photon number are written as 

pi= f P{a')\a'){a'\d^a' (90) 

and 

7o= f P{a')\afd^a', (91) 

respectively, where a' is a complex random variable and P{a') is the classical probability 
density. Consequently, the formulas of the average photon number of amplified signal and 
generated idler beams, Ig and Jj, can be obtained by replacing the term |ap in Eq. fHOj) 
with Iq. 

We then calculate the quantum noise ^IgU) i^ ^'i- p3|) and AJf in Eq. ( l49l) by taking the 
quantum average over the density operator pj. After some algebra, we get the normalized 
intensity noise of the amplified signal (generated idler) beam 

R'si^) = Rsi^)+lsi.){{\aT)-{{\af)r), (92) 

and the normalized intensity difference noise of twin beams 

R't = Rt + %^^ iilaT) - ((l«?))^) , (93) 

h + rli 
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where (la'l") = /P(a')l'^'r'^^'^'5 ^^^d the inequality {\a'\'^) — {{\a'\'^))'^ > holds for the input 
signal with excess noise. In Eqs. f l92|) and fl93|) . the first terms Rs(i) and Rt in the right hand 
sides are the corresponding noise for input signal in a pure coherent state with Jj„ = Jq (see 
Eqs. ( H71) and (1531)): while the second terms in the right hand sides are originated from the 
excess noise of input signal. 

Equations ( p2|) and ( 193|) indicate that the side effect of the excess noise can be eliminated 
under the condition r = Is/U, and the experiment in Ref. J20| has verified this for the case of 
r = 1. In the low gain regime, the value of r = Is/ It is different from ropt for minimizing Rt 
(see Fig. EJ^b)), so the excess noise the input signal is deleterious for the noise reduction of 
Rt. However, in the high gain limit, we have r = Is/h — )■ 1. In this case, eliminating the side 
effect of excess signal noise and minimizing Rt can be achieved simultaneously. Therefore, 
a high gain FOPA is desirable for improving the noise reduction of twin beams. 

VI. SUMMARY AND DISCUSSIONS: 



In conclusion. We have developed the multi-mode quantum theory for analyzing the noise 
characteristics of a non-degenerate phase insensitive FOPA pumped by a pulsed laser. The 
calculation shows that both the noise figure and the intensity difference noise of twin beams 
depend on the JSF. In the high gain regime, the noise figure is generally greater than the 
3dB quantum limit unless the JSF is factorable and the spectrum of the injected signal 
well matches the gain spectrum in signal band; while the intensity difference noise can be 
significantly less than the SNL. To closely resemble the experiments, the quantum noise 
of twin beams generated by a broadband FOPA is numerically studied by taking the real 
experimental conditions into account. In addition to the influences of the quantum efficiency 
of detectors and the excesses noise of input signal, the influence of the collection efficiency 
of the twin beams, determined by the ratio of the pump bandwidth to signal bandwidth p 
and the ratio of the signal bandwidth to filter bandwidth s, are carefully analyzed as well. 
Our theoretical investigation is of practical importance, it not only serves as a guideline for 
optimizing the measured noise reduction of the twin beams, but also helps for understanding 
the quantum noise preformance of pulse pumped FOPA. 

Instead of using the Bloch-Messiah reduction or Schmidt decomposition of the JSF, which 
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has been exploited to theoretically analyze the pulsed CV nonclassical light via the high gain 



spontaneous parametric down conversions in Refs. 23|, |2J], our transformation functions of 
the Bogoliubov transformation are expressed in the infinite series of the gain coefficient G. 
The accuracy of our numerical simulation can be controlled by the two parameters: one is 
the precision of the JSF determined by the step size of the frequency range, the other is the 
order of the series determined by the gain of FWM. So our method has a greater flexibility. 
We believe our theory can be extended to a more generalized OPA, including the degen- 
erate OPA and the phase-sensitive OPA. Moreover, our theory can also be used to study 
the quantum noise of other quantities, such as the noise correlation between the quadrature 
amplitudes of the signal and idler beams. However, it is worth pointing out that the FOPA 
in our model is simplified as a linear amplifier, and we only considered the FWM process. 
Therefore, our calculation only quantitatively explains the results obtained under the con- 
dition of negligible pump depletion [20]. To model the FOPA system more accurately, other 
effects like Raman effect, higher-order FWM and gain saturation effect should be included 
as well 120 1 . 
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